Temperature Dependence of the Effective Bag Constant and the Radius of a Nucleon 

in the Global Color Symmetry Model of QCD 
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We study the temperature dependence of the effective bag constant, the mass, and the radius 
of a nucleon in the formalism of the simple global color symmetry model in the Dyson- Schwinger 
equation approach of QCD with a Gaussian-type effective gluon propagator. We obtain that, as 
the temperature is lower than a critical value, the effective bag constant and the mass decrease 
and the radius increases with the temperature increasing. As the critical temperature is reached, 
the effective bag constant and the mass vanish and the radius tends to infinity. At the same time, 
the chiral quark condensate disappears. These phenomena indicate that the deconfinement and the 
chiral symmetry restoration phase transitions can take place at high temperature. The dependence 
of the critical temperature on the interaction strength parameter in the effective gluon propagator 
of the approach is given. 
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I. INTRODUCTION 



The phase transitions of quantum chromodynamics 
(QCD), for example the evolution between chiral sym- 
metry breaking and its restoration, the color (or sim- 
ply quark ) confinement and deconfinement, have been 
the most active topics in nuclear and particle physics in 
recent years Even though recent investigation has 
provided hints that the phase transitions can be driven 
by the intrinsic characteristics, such as the running cou- 
pling strength and the current quark mass, of the system 
(see, for example, Refs. 04^] ), the more promising and 
much better believed is that the QCD may undergo phase 
transitions into a chirally symmetric and color deconfincd 
phase at high temperature and/or density @- 

To demonstrate the phase transitions, one usually im- 
plements the variation behaviors of not only the fea- 
tures of QCD vacuum and the strong interaction mat- 
ter but also the properties of hadrons at finite temper- 
ature and/or density. On theoretical side, one needs 
in principle QCD, which has been widely accepted as 
the fundamental theory of strong interaction, to carry 
out the investigation. However, as a basic theory, QCD 
still suffers from difficulties in the low energy region, 
which relates directly to strong interaction matter and 
hadrons. Then, besides the approaches of lattice simu- 
lations, QCD sum rules, instanton model(s) and Dyson- 
Schwingcr equations and several models, such as the bag 
model Q, quark- meson coupling model (QMC) [jj-fTTj]. 
Nambu-Jona-Lasinio (NJL) model [l2j, Polyakov-loop 
improved NJL model |l3[, global color symmetry model 
(GCM) [3, and so on, have been developed. The NJL 
model has been widely used since it preserves the fea- 
ture of chiral symmetry and its dynamical breaking and 
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is easy to carry out numerical calculation. Even though, 
with the Polyakov-loop improvement, the quark confine- 
ment effect is included at statistical level, the commonly 
accepted one still only takes into account the point (con- 
tact) interactions among quarks. The bag model is the 
one which handles hadrons as bubbles of perturbative 
vacuum immersed in the physical vacuum. However, all 
nonperturbativc physics is included in a quantity — bag 
constant, which is dealt with a phenomcnological param- 
eter in the model. And the QMC model involves the sim- 
ilar problem. For the GCM, since it can take the result 
of the Dyson- Schwinger equation (DSE) [l5| approach as 
input, it manifests well the properties of chiral symmetry 
and its dynamical breaking. Because the bag constant in 
the model is taken as the difference between the energy 
densities of the perturbative and the physical vacuums, 
the color confinement effect is also handled well in some 
sense. The GCM is then believed to be a quite sophis- 
ticated model which involves as many characteristics of 
QCD as possible. And the NJL model, the QMC model, 
and the bag model can be treated as the special cases of 
the GCM. 

Due to its solid QCD foundation, the GCM has been 
widely taken to study not only the properties of nucleon 
and some mesons [lj, dlHHI hut also the QCD vacuum 
structure [ltl [2(| . It has also been extended to investi- 
gate the properties of strong interaction matter at finite 
temperature and /or density and those of some hadrons 
in the matter |2ll426| . For the property of nucleon in 
finite density matter, it has been discussed with vari- 
ous models for the effective gluon propagator in the DSE 
and the variation behaviors of the mass, the radius and 
the bag constant of the nucleon have been given explic- 
itly j23j . However, in the case of finite temperature, only 
the changing feature of the bag constant has been dis- 
cussed with the Munczek-Nemirovsky model [27[ for the 
effective gluon propagator of the DSE [22j ■ We will then, 
in this paper, discuss some of the properties of a nucleon 
at finite temperature with a sophisticated effective gluon 
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propagator in the DSE. 

The paper is organized as follows. In Sec. II, we de- 
scribe briefly the formalism of the GCM soliton model. 
In Sec. Ill, we describe the algorithm to carry out the 
numerical calculation of the GCM soliton at finite tem- 
perature and the obtained results. Finally we give a brief 
summary and some remarks in Sec. IV. 



II. BRIEF DESCRIPTION OF GCM 

The original action in the global color symmetry model 
(GCM) defined in Euclidean metric is expressed as [fij 



SGCM(q,q) 



d xq(x)(i-f ■ p + m )q(x) 
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+K-ll<f X( fyjZ(x)DZUx - y)jt(y) ,(1) 



where j^(x) = q(x)^-j tl q(x) is the local quark current, 
D ab v {x — y) is the full gluon propagator, m is the cur- 
rent quark mass, g is the quark-gluon coupling constant. 
The Euclidean metric is such that a ■ b = a^b^, and 
{7/i, 7«/} = 2<5 M i,. Taking the gluon propagator to be color 
diagonal in the Feynman-like gauge, i.e., 
8 ab 5 lJLV D{x — y), and applying the Fierz transformation 
to reorder the quark fields, one can rewrite the action as 



and pseudoscalar terms labeled by 9. B e (x, y) are bilocal 
Bose fields. Theoretically, it can be proved that the GCM 
is valid in any gauge even though one takes the Feynman- 
like gauge in deriving the above expression [28j . 
By integrating the quark fields, one gets the action 



S GC M[B B {x,y))] = —Tr \^G~ X \B 6 (x,y)\ 

^B e {x,y)B^{y,x) 



d xd y- 



2g 2 D(x - y) 



, (3) 



where the inverse of the quark propagator can be written 

as 

G- 1 (x, y) = {i 1 -p + m )S(x - y) + A e B e (x, y) . (4) 
Generally, the bilocal fields can be expanded as 



B e (x,y) = B e (x,y)+J2r 6 o(x,y)$ 



(5) 



where the first term is the translation invariant vacuum 
configuration. The second term stands for the fluctua- 
tions of the vacuum which can be identified as effective 
meson fields since the 9 stands for the quantum num- 
ber of Bose fields. In the lowest order, one takes the 
Goldstonc mode, (f)® = {cr, ir}, which is thought of as the 
most important low energy degree of freedom. The vac- 
uum configuration can be determined by the saddle point 



Sgcm[B {x,y))] = d xd yq(x)[(ij ■ p + m )S(x - y) condition 6S G cm/6B$ = 0. One has then 



+A e B 6 (x,y)]q(x) 
dXdV 2g 2 D{x-y) 



B 6 (x,y) = g 2 D(x - y)tr[G( yi x)A e 



(G) 



(2) 



where {A e } are direct products of Lorentz, flavor and 
color matrices of quarks which produce the scalar, vector, 



and the quark self-energy E(x — y) = A 6 B^x^y). The 
equation of quark self-energy in momentum space coin- 
cides with that of the truncated Dyson-Schwingcr equa- 
tion (DSE) 



J 



E(p)= / d A xA e B e {x,y)e^ = g 2 



d 4 q N A a 



I X a 
»7- (q + m) + S( 9 ) 7ai Y 



(7) 



r 



where i M „ = <5 M „ — k^k v /k 2 , with k = p — q, 7 M is the color 
SU(3) matrix. Generally, the quark self-energy function 
can be decomposed as 



Z(p) = S- 1 (j>)-So\p), 
S- 1 {p) = i 1 -p A(p 2 ) + B(p 2 ), 
S o 1 (p) = i l -p + m Q , 



(8) 



Recalling the configuration of the bilocal fields in 
Eq. (5) and considering the Bethe-Salpeter amplitude of 
the mesons and the partial conservation of axial- vector 
current, one can prove [2!| that, when considering the 
most important low energy degree of freedom, i.e., the 
Goldstone mode c/>q = {cr, n}, Eq. (5) can be rewritten as 



and A(p 2 ), B(p 2 ) are scalar functions of p 2 . 
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A e [B e (x,y)-B e (x,y)] = B(x-y) 



(9) 



r 



where _B is just the scalar part of the inverse of the 
quark propagator which can be determined by solving 
the quark's DSE. 

With a nontopological-soliton ansatz[3(i|. the action of 
the GCM soliton with quarks in chiral limit (m = 0) 
can be given fll [i~7l [l8j as 

Sgcm = q{il ■ P ~ a[a(x) + m(x) ■ T7 5 ]}<7 

+ / [f (^) 2 + ^{d^f-Via^zm 



with 
V(a,ir) 



-12 



d 4 P 



In 



A 2 {p)p 2 + (a 2 + Tf 2 )B 2 {p) 



A 2 (p)p 2 +B 2 (p) 



(a 2 +n 2 -l)B 2 (p) 



A 2 {p)p 2 +B 2 (p) 
and the quark meson coupling constant a is given as 

d A p 



(11) 



a(x) 



(27T)" 



-B{p)e 



It is evident that such a quark meson coupling constant 
is just the vacuum configuration of the bilocal fields and 
is independent of the meson fields. 

With the stationary condition of the soliton, one has 
the equations of motion for the quarks and mesons as 



{ij ■ p - a[a(x) + iir(x) ■ T75]} q = , 



V 2 a(f) 



SV 
Sa(r) 



+ Q*(r) = 0, 



(12) 



(13) 



- V 2 7?(r) + + Q* (?) = , (14) 

where Q a and are the source terms contributed from 
the valence quarks, and can be written as 



3=1 



(15) 



Qt((R) = y^V - I d 3 xd 3 yuj(x)B(x - yji^r 

— 1 J J 



3=1 

5(^-R) Uj (y), 
with Zj being the renormalization constant [l] 



(16) 



^(Pl^l^uj® • (17) 



The quark field and er, 7r meson fields can be deter- 
mined by solving Eqs. (12)-(14) self-consistently. As a 
consequence, the corresponding eigenenergies can be ob- 
tained. It is apparent that the meson fields correspond- 
ing to the vacuum configuration can be simply taken as 
a = 1, 7r = due to the (normalized with f 2 ) restriction 
7r 2 + cr 2 = 1. The vacuum configuration is a minimum of 
V(<7,7r) and V(1,0) = 0. In light of the nontopological- 
soliton ansiitz [14J, l3Cf , one can approximate the soliton as 
a chiral bag with bag constant 



V(<r,n)-V(l,0) = 12 



d 4 p 
(2^ 



In 



A 2 (p)p 2 + B 2 (p) 



A 2 {p)p 2 + (a 2 + 7f 2 )B 2 {p) 
I 



(a 2 + n 2 -l)B 2 (p) 
A 2 (p)p 2 + B 2 (p) 



(18) 



With the correction from the motion of center of mass, 
the zero-point effect, and the color-electronic and color- 
magnetic interactions being taken into account, the total 
energy of a bag (involving three valence quarks) is given 
as 



E B (R) = 3sj(R) + -irR 3 



Zo 
R 



(19) 



where Sj (i?) is the energy eigenvalue of the quark's equa- 
tion of motion, R is the radius of the bag, Zq/R denotes 



the contribution of the corrections of the motion of the 
center of mass, zero-point energy, and other effects with 
Zq being a parameter. Just the same as that in Ref. [lij], 
the bag is identified as a nucleon in the present work. 
With the equilibrium condition 

dE B (R) 
dR 

we can obtain the radius of a nucleon. 
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III. ALGORITHM AND NUMERICAL RESULTS 
A. Algorithm 

From the description in last section, we know that the 
property of a nucleon (for instance, its mass, radius and 
bag constant) is determined by the solutions of the equa- 



tions of motion of the quarks and (chiral) mesons in the 
soliton. To solve the equations of motion, one needs 
the solutions A(p 2 ) and B{p 2 ) of the quark's Dyson- 
Schwinger equation. Then after solving the quark's DSE 
in Eq. (7), or more explicitly [with the help of the de- 
composition in Eq. (8)] the coupled equations 



/ 4Si D (P - l)tr [*7 • P t^S(q)T u (p, q)] , 
/ 0^D(p-q)tr[t^S(q)T„(p,q)], 



(20) 



where Cf is the eigenvalue of the quadratic Casimir oper- 
ator in the fundamental representation of the color sym- 
metry group [for SU(N C ), C F = (N 2 - l)/2N c , it reads 
4/3 at N c = 3]. 

To investigate the temperature dependence of the 
property of a nucleon, one should at first discuss the form 
of the quark's DSE at finite temperature. 

It has been well known that the appearance of 
(nonzero) temperature T in the QCD reduces the 0(4) 
symmetry to 0(3). Thus the quark's four-momentum 
p should be rewritten as p — (p,ui n ), where w„ = 
(2n + 1)tvT (n G Z) are the discrete Matsubara frequen- 
cies of the quark, and the four-dimensional integral needs 
to be replaced by [3l[ 



d 4 P , rp f d 3 p 



(27T) 



(2tt)= 



(21) 



The decomposition of the dressed quark propagator needs 
to be rewritten as 

S~\p,u,) = i ; ?j?A(\p\,u) r ) + i'y 4 (j n C(\p\, ujJ+B(\p\,ujri) ■ 

(22) 



Furthermore, the gluon propa gato r at finite temperature 
can be generally expressed as |31| 



with the transverse and longitudinal projectors 



- J (k 



0. 



fi or/and v = 4, 



(23) 



(24) 



-Pl(k). 



Due to the lack of detailed information about the gluon 
propagator at finite temperature, one usually naively as- 
sumes that the transverse and longitudinal parts of the 
gluon propagator are equal and independent of temper- 
ature (see, for example, Ref. |32|), i.e., one has approxi- 
mately DT{k,uj n ) = Dj J (k,cj n ) = D(k). In this assump- 
tion and bare vertex approximation, we get three coupled 
integral equations for the functions A(p, w n ), C(p,u> n ) 
and B(p, ui n ) as 



A(\p\,cj n ) = l + g 2 T%Y. I T^H^X {HP - « )fc2 + ' k)]A{\q\,u m ) + 2w m £l k {p ■ k)C{\q\,oj m )} , 

C(\p\,uj n ) = l + g 2 T^J2 f 1" {^ nUJmk2 + 2uj n uj m nl}C(\q\,io m ) + 2w n n k (p-k)A(\q\,u m )} , (25) 

n m 

B(\p\,u Jn ) = mo + g 2 TC F J2 I &^- W (W> u ™)> 

i 

where k 2 = k 2 + f2fe = oj n — LJ m , and A = q 2 A 2 + As mentioned in last section, in view of the non- 

ijJmC 2 + B 2 . topological soliton ansatz, one can take a nucleon as a 
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soliton bag. In the chiral limit (mo = 0), one has found 
that there exist two types of solutions for the quark's 
DSE. One is the Nanmbu-Goldstone solution which cor- 
responds to the chiral symmetry spontaneously broken 
phase. The other is the Wigner solution which represents 



the state with the chiral symmetry One gets then the 
effective bag constant as the pressure difference between 
the Nambu-Goldstone solution and the Wigner solution, 
which reads 



3{T) = P[Q 



NG] 



4^E 



P[Q 

d 3 p 
(2^ 



In 



A 



NG 



Mr 



j^Ang + uSLC NG p*A w + uJLCw 



A 



JVG 



Ai 



(26) 



where A NG = fA 2 NG + io 2 n C 2 NG + B 2 NG , A w = fA 2 w + 
^0%, and A NG , C NG , B NG , A Wl C w denotes the 
Nambu-Goldstone solution, the Wigner solution for the 
DSE, respectively. 



With the solutions of the DSE as input, we can have 
the explicit expression of the quark's equation of motion 
at finite temperature as 



[ij-pA(f, T) + i l4 u n C{p, T) + B(p,T)] Uj (p, T)- 



d 3 k 



-B 



,T) a{p-k) + i l5 T-Tf(p-k) Uj (k,T) = 0, (27) 



r 



where a = a — 1 . After solving the related eigencquation 
we can obtain the eigenenergy of the quark at finite tem- 
perature, Ej (T) . One can in turn study the properties of 
a nucleon at finite temperature by extending the GCM 
soliton model described in last section. 

In the case of zero temperature, one usually takes only 
the lowest energy for the £j(R) in Eq. (19). At finite- 
temperature, due to the influence of temperature, the 
contribution of the excited states of the quarks must be 
included. Then the energy of the bag, i.e., the mass of a 
nucleon, at finite temperature T should be written as 



M(T)=E B (T) = 3e j (T) 



Zo 
R 



--kR 3 SS{T) , (28) 



where £j{T) is the average of a quark's energies at all 
possible states. 

In the spirit of the most simple approximation of the 
GCM (3] (mentioned at the end of the last section), 
the quarks in the bag (GCM soliton) at finite tempera- 
ture can also be regarded as the free one which satisfies 
the Dirac equation, and the energy eigenvalue can be ex- 
pressed as 



,-cn 



R{T) ' 



(29) 



where j denotes the quantum number labeling the energy 
level. Since quarks are fermions, we take the Fermi-Dirac 
statistics to evaluate the average energy of each quark in 



the bag at finite temperature, and have 

+ e e 3 (T)/T ' 



(30) 



3=0 



where N is the degeneracy of quarks. 



Then, by solving the stability condition dM ^ 



dR 



i.e. 



dE B {T) _d£j(T) Zy 



dR 



dR 



+ AitR'33{T) = , (31) 
R* 



we can obtain the (stable) nucleon radius R{T) and the 
mass of a nucleon (the energy of the bag) M(T). 

It is apparent that after solving the quark's DSE and 
in turn the equations of motion of the quark and meson 
fields at zero and nonzero temperature, we can obtain 
the bag constant, the mass and the radius of a nucleon 
in the corresponding circumstance and discuss the vari- 
ation characteristic of the property with respect to tem- 
perature. To solve the quark's DSE, we take a simplified 
form of the effective gluon propagator in Ref. [33| 



(32) 



g 2 D(k) = A-n 2 D—s exp =■ 



where D and u> are dimensional parameters that can 
be determined by fitting empirical date. Such an ef- 
fective gluon propagator is naturally an extension of 
the Munczck-Nemirovsky model [27J and consistent with 
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those given in lattice QCD calculations (see, for instance, 
Refs. |34l . HU) and solving the coupled DSEs of quark, 
gluon, and ghost (see, for instance, Ref. [H]). It has also 
been shown to be successful in describing many hadron 
properties [H, [H, UMl • 



B. Calculation and Numerical Results 

We first solve the quark's DSE at chiral limit (m = 0) 
and zero temperature with parameters D = 1.0 GcV 
and lo = 0.5 GeV, which have been widely used (see, 
for example, Refs. [H, Ull)- The obtained result of the 
Nambu-Goldstonc solution of the DSE at zero tempera- 
ture is displayed in Fig. 1. It shows evidently that our 
result reproduces exactly that given in Ref. [371 ] . 




10' 

p 2 [GeV 2 ] 



FIG. 1: (Color online) Calculated result of the Nambu- 
Goldstone solution of the quark's DSE at zero tempera- 
ture, with parameters in the effective gluon propagator u — 
0.5 GeV, D = lGeV 2 . 

We then solve the quark's DSE at nonzero temper- 
ature with the same effective gluon propagator. Fig- 
ure 2 illustrates the obtained results of the functions 
A(\p\,T), C(\p\,T), B(\p\,T) in the Nambu-Goldstonc 
solution and the mass function M(\p\,T) at a tempera- 
ture T = 30 MeV, as an example of those at nonzero tem- 
perature. We find from the figure that, just as expected, 
as the temperature is lower, the functions A(\p\,T), 
C(\p\,T), B(\p\,T) and M(\p\,T) have the correct zero 
temperature limit. Hence, even though we have not in- 
cluded explicitly the temperature effect in the effective 
gluon propagator, the calculated result can demonstrate 
the temperature dependence of quark propagator with 
quite high precision. In addition, it can be noticed from 
Fig. 2 that functions j4(|p|,wo) and C(|p|,wo) have the 
same behavior when the temperature is low, just as that 
of the quark propagator at low chemical potential [40j |. 
To study the Matsubara frequency dependence of the 
quark propagator, we illustrate the variation behaviors 
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FIG. 2: (Color online) Calculated results of the functions 
A(\p\,T), C(\p\,T), B(\p\,T) and M(\p\,T) in the Nambu- 
Goldstone solution of the quark propagator at a temperature 
T = 30 MeV, with parameters in the effective gluon propaga- 
tor u = 0.5 GeV, D = 1 GeV 2 . 



of the functions A(|p|,cj„) and B{\p\,uj n ) at a tempera- 
ture T — 30 MeV in Fig. 3 [since Fig. 2 has shown that 
the functions A(\p\,u n ) and C(\p[,u n ) at T = 30 MeV 
are almost exactly equal to each other, we do not show 
the function C(|p|, as a representative. From Fig. 3, 
we can find that, when the temperature is low, both func- 
tions A and B involve an approximate bilateral symme- 
try about n [since for a fixed temperature, u„ is propor- 
tional to n due to the definition u n = (2n + l)7rT], and 
the function B decreases rapidly as the Matsubara fre- 
quency increases. It gives us a posterior knowledge that 
when we carry out the summation of the Matsubara fre- 
quencies w„, it is not necessary to do that up to a very 
large number of n. 

To demonstrate the temperature dependence of the 
quark propagators explicitly, we illustrate the calculated 
variation behaviors of the Wigncr solution at zero mo- 
mentum and zero mode Matsubara frequency [i.e., the 
functions Avp([pj = 0,cjo) and CV([p| = 0,wo)] and the 
Nambu-Goldstone solution under the same conditions 
[i.e., the functions Ang{\p\ = 0,uq), Cng{\p\=0,uo) and 
Bng{\P\ =0,wo)] with respect to temperature in Fig. 4 
as a representative. From Fig. 4, we can find that, when 
the temperature is low, the functions Amg and Cng co- 
incide with each other very well. So do the Aw and 
Cw except that the temperature for the deviation be- 
tween them to appear is lower. It indicates that only as 
the temperature is quite high, the breaking from 0(4) 
symmetry to 0(3) symmetry (especially, for the physical 
state, i.e., the Nambu-Goldstone state) becomes obvious. 
Moreover, the decreasing feature of the function -Bjvg is a 
manifestation of gradual restoration of chiral symmetry. 

With the solutions of the quark's DSE as input, we 
solve the equations of motion of the quark and the chiral 
meson fields, and then obtain the property of a nucleon 
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2010 IPltGeV] 




10 10 20 

n 



FIG. 3: (Color online) Calculated results of the functions 
A(]p|,aj n ) and B(\p\,u} n ), where cj„ = (2n + 1)ttT, in the 
Nambu-Goldstone solution of the quark propagator in chiral 
limit and at a temperature T — 30 MeV [(a), (b), respectively] 
and the special case at |p| = [(c), (d), respectively]. The cal- 
culations are also carried out with parameters u> = 0.5 GeV, 
D = 1 GeV 2 in the effective gluon propagator. 



at zero and nonzero temperature. The obtained prop- 
erty of a nucleon at zero temperature is bag constant 
^(0) = (162 MeV) 4 , mass M(0) = 939 MeV (with the 
parameter Zq is fixed as 3.08), and radius R(0) = 0.85 fm. 
It should be noted that the presently fixed value of the 
parameter Zq, 3.08, is larger than the usually taken one, 
1.84. Such a large value arises from the fact that the 
term —Z /R is a combination of the contributions from 
not only the zero-point energy but also those of the color- 
electronic and color-magnetic interactions, the correction 
on the motion of center-of-mass and other effects. And it 
is consistent with the most recent result [4l| and our pre- 
vious results [23|. The gained variation behaviors of the 
nucleon's bag constant, mass, and radius with respect 
to temperature are illustrated in Figs. 5, 6, 7, respec- 
tively. From Figs. 5-7, one can notice that, with the 
increasing of temperature if it is below a critical one, 
the bag constant and the mass of a nucleon decrease, 
and the radius of the nucleon increases. At critical tem- 
perature T = 133 MeV, the bag constant and the mass 
of the nucleon vanish and the radius tends to be infi- 
nite. It manifests that the nucleon can no longer exist 
as a bag soliton, so that the quark deconfincmcnt hap- 
pens. Admittedly, such a obtained critical temperature 
may be model-dependent, however, the gradual variation 
features of the nucleon's property indicate that the de- 
confinement phase transition process is that, with the in- 
crease of temperature, the nucleons in the matter touches 
with each other at first due to the increase of the radius, 
then the fields of the ingredients of the nucleons mixed 
with each other and the bound strength gets weaker si- 
multaneously. As the bound (the bag constant) vanishes, 
the deconfincmcnt phase transition occurs. Therefore, 
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FIG. 4: (Color online) Calculated results of the temperature 
dependence of the Wigner solution of the quark's DSE at 
zero momentum and zero mode Matsubara frequency [i.e., 
the functions Aw(|p| = 0,a>o) and CV(|pl = 0, Wo)] (panel 
(a)) and the Nambu-Goldstone solution under the same con- 
ditions [i.e., the functions Ajvg(|p| = 0,cj o )i Cng(\p\ = 0,cj o ) 
and Bng(\p\ — 0, lj )] (panel (b)). The calculations are also 
carried out with parameters u> — 0.5 GeV, D — 1 GeV 2 in the 
effective gluon propagator. 
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FIG. 5: Calculated result of the variation behavior of the bag 
constant of a nucleon with respect to temperature. 



the temperature driven deconfinement process may be, 
in fact, a crossover but not a low order phase transition. 

When discussing the QCD phase transition, one usu- 
ally interests in the chiral symmetry and its dynamical 
breaking, too, and takes the chiral quark condensate as 
a order parameter in the case of chiral limit, which is 
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FIG. 6: Calculated result of the variation behavior of the mass 
of a nucleon with respect to temperature. 
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FIG. 8: Calculated result of the variation behavior of the 
chiral quark condensate with respect to temperature. 
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FIG. 7: Calculated result of the variation behavior of the 
radius of a nucleon with respect to temperature. 

defined as 

n =+°° r ^3 

-(qq)=N c TY / / 7A3 tr[S(p, 0J n )} , (33) 

We then calculate the temperature dependence of the 
chiral quark condensate. The obtained result is shown 
in Fig. 8. The figure displays evidently that the chiral 
quark condensate decreases gradually with the increase 
of temperature and vanishes at a critical temperature. It 
indicates that dynamical chiral symmetry breaking effect 
gets weaker and weaker with the increase of temperature 
and the chiral symmetry can be restored as the tempera- 
ture reaches the critical one. In the case with parameters 
lo = 0.5 GeV, D = 1 GeV 2 in the effective gluon propa- 
gator, the critical temperature for the chiral symmetry 
to be restored is also approximately 133 MeV, the same 
as that for the deconfinement. 




D' 



FIG. 9: Calculated result of the variation behavior of the 
scaled critical temperature with respect to the scaled strength 
parameter D' in the effective gluon propagator, where T' c = 
T c /u and D' = D/lj 2 are the scaled quantities. 



As mentioned above, the critical temperature for the 
QCD phase transitions to happen may be parameter de- 
pendent. To demonstrate the parameter dependence ex- 
plicitly, we carry out a series calculations with various 
values of the coupling strength parameter D and the 
screening width parameter lo in the effective gluon prop- 
agator. Due to the good behavior of the Gaussian-type 
gluon propagator, we can scale the T c and D by to with 
definition T' c = T c /u and D 1 = D/uo 2 . The obtained 
variation behavior of the T' c with respect to the D' is 
displayed in Fig. 9. One can find easily from the figure 
that the critical temperature increases when the coupling 
strength gets larger. Furthermore, there exists a critical 
scaled coupling strength D' , below which the critical tem- 
perature for the deconfinement maintains zero. In fact, 
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below the critical coupling strength, the quark's DS equa- 
tion does not have a Nambu-Goldstone solution [H, Q . In 
other words, there exists only quarks in chiral symmetry. 
In turn, we have only deconfined quarks, but no nucleons 
(more general, hadrons) even if the temperature is zero. 

IV. SUMMARY AND REMARKS 

In this paper we have calculated the temperature de- 
pendence of the quark propagator by solving the quark 
DSE with a Gaussian-type effective gluon propagator. 
Based on the calculations, we investigated the tempera- 
ture dependence of the bag constant, the mass and the 
radius of a nucleon in the framework of the GCM soli- 
ton model. It shows that, as the temperature is lower 
than a critical value, the bag constant and the mass de- 
crease and the radius increases with the increasing of the 
temperature. In the case with parameters ui = 0.5 GeV, 
D = 1 GeV 2 in the effective gluon propagator, the critical 
temperature is found to be about 133 McV. At the criti- 
cal temperature, the bag constant and the mass decrease 
to zero and the radius increases to infinity. It means that 
the nucleon can no longer exist as a bag soliton, so that 
the deconfinemcnt phase happens. This indicates evi- 
dently that the quark deconfincment phase transition can 
take place at high temperature. Moreover, we give the 
dependence of the critical temperature on the interaction 
strength parameter D in the effective gluon propagator. 
It shows that, as the interaction strength parameter is 
larger than a critical value, the critical temperature in- 
creases with the increasing of the strength parameter. 

Even though the temperature dependence of some of 
the properties of nucleon is given with some approxi- 
mations and model parameters in the present work, the 
qualitative behavior would be universal and it is the first 
one given with quite a sophisticated approximation of 
QCD. Of course, there are various aspects to be im- 
proved. For example, we take the commonly used ef- 
fective gluon propagator, which is independent of tem- 
perature, to solve the quark DSE and make use of the 



preliminary GCM soliton model [l4|, HH, [H| . It is neces- 
sary to implement the real GCM soliton model [13, [H| 
with solving at first the coupled DSEs of the quark, gluon 
and ghost, and then the coupled equations of the quark 
and the chiral fields, with the inclusion of the temper- 
ature effect. In more detail, for the quark gluon inter- 
action vertex, we take simply the bare vertex 7 M in our 
present work, in fact more realistic vertex functions, such 
as the BC vertex [42| , which has been shown to be able to 
improve the calculation of meson properties greatly [3!| , 
even the BC vertex together with the transverse part 
being included simultaneously [43T4451 ]. should be imple- 
mented to make the calculation with much more solid 
QCD foundation. It is also necessary to notice the dimen- 
sionless quantity Zq which takes account of the contri- 
butions of the zero-point effect, the color-electronic and 
color-magnetic interactions, the motion of center-of-mass 
and all the others. In our present work, we handle it, in 
the commonly taken way, as a free parameter to be fixed 
by the property of a nucleon in free space. In fact, all the 
aspects of the Zq are quite complicated and have recently 
been paid great attentions (see for example, Refs. (46l.l471| 
tried to evaluate the zero-point effect part from the gluon 
field fluctuations directly). Furthermore, extending the 
result of the thermal Casimir effect in ideal metal rect- 
angular boxes HH, we infer that the Zq (at least, the 
zero-point effect part) may depend on temperature. It 
would then be interesting to study the temperature de- 
pendence of the parameter Zq . The related investigations 
are in progress. 
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